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Abstract 

We introduce the Callan-Symanzik method in the description of anisotropic as well as isotropic 
Lifshitz critical behaviors. Renormalized perturbation theories are defined by normalization con- 
ditions with nonvanishing masses and at zero external momenta. The orthogonal approximation 
is employed to obtain the critical indices t]l2, Vl4 an d vla diagramatically at least up to 

two-loop order in the anisotropic criticalities. This approximation is also utilized to compute the 
exponents vla and in the isotropic case. Furthermore, we compute those exponents exactly for 
the isotropic behaviors at the same loop order. The results obtained for all exponents are in perfect 
agreement with those previously derived in the massless theories renormalized at nonzero external 
momenta. 
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I. INTRODUCTION 



Renormalized perturbative expansion of a X(p A field theory is the natural mathematical 
setting in statistical mechanics to determine the critical properties taking place in second or- 
der phase transitions. For ordinary critical systems, the theory is defined on a (^-dimensional 
Euclidean space and has quadratic kinetic terms in momenta space, as usual in quantum 
field theory. These criticalities can be generalized by adding certain combinations of m quar- 
tic derivatives of the field to the kinetic term of the bare Lagrangian density. In fact, when 
the m coefficients of the usual quadratic derivatives of the field vanish, the Lagrangian den- 
sity describes the so-called m-axial Lifshitz critical behaviorsQ, [j, [3]. Physically, the bare 
Lagrangian describes an Ising model with ferromagnetic exchange forces between nearest 
neighbor spins as well as antiferromagnetic interactions between second neighbors along m 
space directions. The competition among ferro- and antiferromagnetic couplings of the spins 
in the model induces the vanishing of the m terms in the quadratic derivatives of the field 
keeping, however, the quartic derivatives along the competing directions. The m higher 
derivative terms will generate space anisotropy whenever m < d, with two inequivalent 
subspaces, whereas if m = d there is only one isotropic space. 

From the quantum field-theoretic side, interest in the model stems from the fact that it 
has a natural connection with Lorentz symmetry breaking for fields propagating in specific 
backgrounds in the long-distance limit. Indeed, a recent proposal postulates the existence of 
a physical fluid originating from a scalar field whose vacuum expectation value changes with 
a constant velocity [4]. This background breaks the time diffeomorphism symmetry generat- 
ing a type of Goldstone boson for this breaking of Lorentz invariance, the ghost condensate 
(analogous to the Goldstone bosons which arise from spontaneous breaking of internal/gauge 
symmetries). It mimics a type of dark energy, whose low-energy effective Lagrangian in flat 
spacetime has kinetic terms with quadratic time derivatives and quartic space derivatives. 
The ghost condensate mixes with gravity in a kind of Higgs mechanism giving rise to a 
nontrivial modification of gravity in the infrared [5j. The competing interactions between 
attractive and repulsive components of gravity at large scales manifest themselves in the 
absence of quadratic space derivatives of the ghost condensate in the kinetic Lagrangian 
implying that this flat spacetime is not the conventional Minkowski space. After a Wick ro- 
tation in time for this quantum field theory, the kinetic Lagrangian for the ghost condensate 
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(tt) can be identified with that for the order parameter field (magnetization <f>) in the m-axial 
Lifshitz critical behavior for m — 3. Therefore, figuring out the perturbative structure of 
Lifshitz points is worthwhile in order to gain a better comprehension of this kind of Lorentz 
breaking quantum field theory as well as the large distance effects on the gravity sector in 
this model. 

Renormalization group and e-expansion ideas (e = 4 — d) were originally introduced 
to provide the perturbative determination of critical exponents using diagramatic methods 
in momenta space [6(. A reformulation of this method was introduced to compute the 
critical exponents in a massless theory renormalized at nonzero external momenta 3, |8[]. The 
associated renormalization group equation along with further developments in evaluating 
Feynman diagrams have been adapted to study critical exponents of the Lifshitz type when 
more than one characteristic length scale are present Critical exponents of Lifshitz 
points have previously been obtained applying these ideas in massless theories renormalized 



at nonzero external momenta 
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An alternative to the above methods is to formulate the problem in a massive theory. 
This permits to comprehend under what conditions the scale invariance of the solution to 
the renormalization group equation for the renormalized 1PI vertex parts is guaranteed. In 
addition, this formulation is convenient to treat renormalization issues of arbitrary quantum 
field theories and it is a healthy test for the universality hypothesis, which guarantees that 
universal quantities (e.g., the critical exponents) do not depend on the scheme of renor- 
malization. For ordinary critical systems, the asymptotic scale invariance of the theory 
in d = 4 — e was established more rigorously after the first developments Ref. js] through 



the use of Callan-Symanzik equations [12|, [13J in the explicit computation of the critical 
exponents for a massive theory [14]. However, the same approach has not been applied for 
investigating more general critical systems so far. In fact we would like to know if such a 
method could be adapted to critical systems exhibiting more than one characteristic length 
scale, such as the m-fold Lifshitz critical behavior. It is the purpose of the present work to 
show that such a mission can be accomplished, following closely the arguments put forth in 



14] . Here, we derive the critical exponents rj and v in the massive case by introducing an 



appropriate Callan-Symanzik formalism with two independent mass scales in the bare La- 
grangian corresponding to the correlation lengths £t 2 and associated to the anisotropic 
critical behaviors. The massive theories are renormalized at zero external momenta. A 
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similar treatment will be shown to be valid for the isotropic behaviors with only one mass 
scale. Consequently, we derive the critical exponents characterizing the anomalous dimen- 
sions of the field 4> and composite operator (p 2 , respectively, within an e^-expansion at least 
at two-loop order using Feynman's path integral method in momentum space. 

The anisotropic diagrams are calculated using the orthogonal approximation introduced 
in [h]]. The isotropic diagrams are more flexible and can be calculated simultaneously 
using the orthogonal approximation as well as exactly. The exact computation was recently 
demonstrated in the massless theory 15| and shall be calculated here in the massive theory. 
The exponents computed diagrammatically in the present work agree with those obtained 
in the massless theory renormalized at nonvanishing external momenta as expected. 

The organization of the paper is as follows. In Sec. II we present the normalization 
conditions for the anisotropic and isotropic Lifshitz critical behaviors in the massive case. 
The Callan-Symanzik equations and Wilson functions are defined in Sec. III. The critical 
exponents for the anisotropic cases are presented in Sec. IV. The computation of critical 
indices for the isotropic cases is the subject of Sec. V. The conclusions and perspectives of 
the present work are discussed in Sec VI. The relevant Feynman integrals are calculated in 
the appendixes. Appendix A contains the information concerning the computations involved 
in the anisotropic behavior using an approximation to resolve the Feynman integrals. The 
isotropic integrals are computed using the same approximation in Appendix B. Finally, the 
exact isotropic integrals computed exactly are summarized in Appendix C. 



II. NORMALIZATION CONDITIONS FOR THE MASSIVE THEORIES 

Before setting the normalization conditions let us recall some basic facts about the Lifshitz 
critical behaviors. The bare Lagrangian density is that of a self-interacting scalar field given 
byQ : 

L=^\V 2 m<Po\ 2 + \\ V(d-m) 00 I' + S -\ Vm 00 |' + ^00 + JjV0Q- (1) 

The Lifshitz point, where fi T = (for the temperature fixed at the Lifshitz value T = Tl) 
and 5q = is the multicritical point of interest in this framework. Consider the anisotropic 
criticalities. The first term in the above expression should be multiplied by a constant in 
order to make sense on dimensional grounds. Fortunately, the condition 5q = 0, which will 
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be implemented in all the subsequent discussion, permits the dimensional redefinition of 
momentum scales along the competing m-dimensional axes with quartic momentum depen- 
dence. They have half the value of a genuine momentum scale and we can get rid of the 
constant in front of this unusual kinetic term. In the above Lagrangian we introduce two 
independent bare masses in order to perform the appropriate scale transformations con- 
sistent with the two independent correlation lengths characterizing these situations. They 
have different powers in order to emphasize that the momenta subspaces are inequivalent. 
Thus, the bare mass /i2, as well as its renormalized counterpart «2 have half the value of 
a genuine mass scale. These distinctions do bring new insights for it is convenient that 
mass and momentum have the same canonical dimension in each subspace separately. In 
addition, two independent cutoffs are needed to provide independent flows in the parameter 
space as far as the masses are concerned. Together with two coupling constants chacterizing 
each subspace, we have an apparent overcounting which allows the renormalization group 
transformations to be treated independently in the inequivalent subspaces. In the isotropic 
critical behavior there are some differences: the second term in the Lagrangian is absent 
and there is only one type of bare (and renormalized) mass(es). 

It is interesting to consider the renormalized theories that can be constructed out of the 
Lagrangian (1). To do so, we shall follow the conventions addopted in 10] . The reader 
should also consult Ref. [8| for further notation details. 

Let us turn our attention to the anisotropic cases. The two subspaces, namely one with 
a quadratic dependence on the momenta and the other with a quartic dependence on the 
momenta are orthogonal to each other and can be treated independently. In fact, we can 
define two sets of normalization conditions associated to either subspace as follows. The 
quadratic subspace is defined by fj,2,^2 = 0, fii,X\ ^ 0. The latter induce renormalized 
quantities mi and g\. Since no infrared divergence takes place in massive theories, we 
can renormalize the theory at zero external momenta. The one-particle irreducible(lPI) 
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renormalized vertex parts are defined through the following normalization conditions: 

r R(l)(°' m l'5 , l) = m l' ( 2a ) 

dT^l (p, mi , gi ) 

df2 Ip 2 =° = ( 2b ) 

Tf {1) (0,m ugi )= gi , (2c) 

rg' 1 1 ) ) (0,0,m 1 ^ 1 ) = l, (2d) 

rfV(0,m ugi ) = 0. (2e) 

The subscript 1 makes explicit reference to the quadratic or noncompeting subspace. The 
most general vertex parts at this subspace have nonvanishing external momenta along the 
quadratic directions as well as zero external momenta components along the quartic direc- 
tions. 

Similarly, the definition of renormalized theories along the quartic or competing directions 
can be performed by setting first all the external momenta along the noncompeting directions 
equal to zero and taking Ax = with ^2,^2 7^ 0. The resulting 1PI renormalized parts in 
the competing subspace can be rendered finite by using the following set of normalization 
conditions: 

r i?(2)(°' m 2, ^2) = mj, (3a) 
dT ( ^ 2) (p,m 2 ,g 2 ) 



dp 



A 



(3b) 



r 2(2)(°' m2 '# 2 ) = 92, (3c) 
rg$(0,0,m2,(/ 2 ) = l, (3d) 
rg (2 2 ) ) (0,m 2 ,^ 2 ) = 0. (3e) 

The subscript 2 refers to the competing subspace. Therefore, the competing and noncom- 
peting subspaces are effectively separated out. In order to be completely similar to the 
massless case normalization conditions devised in Ref . [lo| . we can fix the mass scale of the 
two-point functions in either subspace by the choices vr?J = 1. Furthermore, the isotropic 
critical behavior can be defined using only the second set of normalization conditions from 
the anisotropic behaviors, with different values of renormalized mass and coupling constants, 
say, 777.3,(73. Next, we apply these conditions to the discussion of the Callan-Symanzik equa- 
tions for the anisotropic and isotropic m-axial critical behaviors. 



III. THE CALLAN-SYMANZIK EQUATIONS 



Let the label r = 1,2 refer to the different external momenta scales involved in the 
general Lifshitz critical behavior, as discussed above for different normalization conditions 
in the anisotropic and isotropic cases. In terms of the bare quantities, the renormalized 1PI 
vertex parts are defined by 

r{ R(r)\Pi(T),Qi(T),9r,m T ) = Zj iT) Z^ {T) (T iN ' L) (p i(T) ,Q i{T) ,X T ,fi T ,A T ) (4) 

-^o^r^'^Q^), Q( T ), K, Vt, K)\q 2 (t) =o) 

where p^ T ) (i = 1, N) are the external momenta associated to the vertex functions F^ffi 
with iV external legs, Qi< T ) {i = 1,...,L) are the external momenta associated to the L 
insertions of (f> 2 operators and A T are cutoffs characterizing each inequivalent subspace. In 
terms of the dimensionless couplings u T , the renormalized and bare coupling constants are 
given, respectively, by g T = u T (m 2r )~^ , and A T = %(mJ T )^, where = 4 + y — d. In 
addition, we can write all the renormalization functions in terms of u T . 

By expanding the dimensionless bare coupling constants u OT and the renormalization 
functions Z^ T ), = Z^Z^^ in terms of the dimensionless renormalized couplings u T 

up to two-loop order as 

u T = u T (l + a\ T u T + a 2T ul), (5a) 
Z^ T ) = 1 + b 2T u 2 T + h T u 3 T , (5b) 
(r ) = 1 + c lT u T + c 2t u 2 t , (5c) 

we can unravel the renormalization structure of the Callan-Symanzik (CS) equations as they 
arise in the two situations. Let us now make a distinction in the two cases by considering 
separately the CS equations for anisotropic and isotropic spaces. 



A. Anisotropic 

In spite of having two independent mass scales, the asymptotic behaviors of the 
anisotropic cases at the critical dimension d c = 4 + y proceed along the same lines as 
the ordinary critical behavior at 4 dimensions described by a A0 4 field theory. Therefore, 
we shall limit ourselves to the situation where d = 4 + y — e L . 
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Before deriving explicitly the equations which relate theories renormalized at different 
mass scales, let us recall some useful definitions for certain quantities in terms of dimen- 
sionless coupling constants. The feeta-functions and Wilson functions for the anisotropic 
behaviors are defined by: 

du 

Pr = (m T -^-), (6a) 

T*(r)(Ur)=/a r ^^ (6b) 
l^{r){U T ) = -P r 7^ (6C) 

are calculated at fixed bare coupling A r . Since the theory is massive, scale invariance will 
be achieved in the limit where the cutoffs A T go to infinity. The /3 T -functions can be cast in 
a more useful form in terms of dimensionless quantities, namely, 

Pr = -re L (^-) . (7) 
By differentiating Eq.(4) with respect to lnm T we find 



t^tT + 0r^~ - 7^7<Kt)(«t) + L^2 {T) (u T ))T { ^ L } \p i{T) ,Q i{Th g T ,m T ) 



[m. 

c)m T c)g T 2 



-$N,oh,2(K T T ) 2 B T = m T T (2r - ^(t))^^ (Pi(r), Qi(r), 0,9t, (8) 

where B T is a constant used to renormalize and 7^2( T )(u T ) = — j3 T — q^ t) ■ The last 

equations are the Callan-Symanzik equations for anisotropic vertex parts with arbitrary 
composite operators. The right hand side can be asymptotically neglected order by order 
in perturbation theory. Consequently, the Callan-Symanzik (CS) equations have asymp- 
totically the same form as the renormalization group (RG) equations previously derived in 



Ref. 



io|. 



For the sake of simplicity, consider the vertex parts without composite operators with 
L = 0. A dimensional redefinition of the momentum components along the competing 
subspace has been performed (see the discussion following Eq.(l) in the last section). Since 
the resulting "effective" space dimension for the anisotropic cases is d e ff = d — y, where d 
is the space dimension of the system, dimensional analysis yields 

Pt 



s 



Now consider the asymptotic part of the corresponding vertex function satisfying the 
homogeneous CS equation. The solution is given by: 

rW (i. m M r N r ipT) < >( \\J< 1 /q\ 

T asR(T)(hT)> u T,—) = exp[-— j ^ {t) (u t ( Pt ))——1 (9) 

r i^(r)(^(r)^r(Pr),m r ), 



where 



J U T 



Ut{Pt) du' 



From now on we shall drop the subscript for asymptotic in the vertex parts satisfying the 
homogeneous Callan-Symanzik equations. The values of the coupling constants u TOO yielding 
the eigenvalue conditions (3 T (u TOO ) = do exist, and precisely at this point the solutions to 
the CS equations under a scale in the external momenta obey the equations 

r R(r){Pr k i(r),UTOo,m T ) = p T (11) 

The dimension of the field is defined by 

CWwI-^^W. (12) 
This in turn implies that the definition of the anomalous dimension of the field through 

d— — 

d 4>(r) = —2*- ~ 1 + 27 corresponds to rj T = 7^( t )(m TO o)- 

If we consider the vertex parts including composite operators, we can identify the anoma- 
lous dimension of the composite operator 2 as follows. For u T = u TOO the aymptotic behavior 
of the vertex parts are given by ((N, L) ^ (0, 2)) 

r RiT ){pTki(T),pTPi(T),u TOO ,m T )=p T (13) 

Writing the coefficient in the right hand side as p}^ 2 ^ ^ 2(t) , we conclude that 

g^> 2 (t) — — 2t + 70 2 ( T )( M roo)- The correlation length exponents can now be obtained through 
the relation v~ x = —d^^ = 2r — 702( t )(m too ). We now have the resources to calculate these 
two exponents diagramatically. 
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B. Isotropic 

Let us briefly discuss the isotropic case d — m. Owing to the existence of only one 
type of scaling transformation, all we have to do is to concentrate on the competing sector 
(r = 2) of the anisotropic treatment performed above. That means that we start with the 
mass scale m-s, dimensionless coupling constant w 3 and subscript 3 in all vertex functions in 
order to avoid confusion with the competing subspace of the anisotropic cases. There are 
3 main differences. First, the critical dimension is 8. Therefore, the expansion parameter 
is e L = 8 — m. Second, the effective space dimension is y. Finally, the beta function does 
not have the global factor of two as in the competition directions of the anisotropic case. 
Instead, in terms of dimensionless parameters one has = -eL(^f)" 1 - Following the 
same trend as before, we conclude that the anomalous dimension of the field in this case 
is i] Li = 773 = 7^(3) (m 3oo ), where w 3oo is the value of the coupling constant which yields 
/5 3 (w 3oo ) = 0. In addition, the anomalous dimension of the composite operator is related to 
the correlation length exponent through the relation v^ 1 = —d^^ = 4 — 74,2(3) (w 3oo ). Let 
us turn our attention to the calculation of these critical exponents by means of diagramatic 
expansions. 

IV. ANISOTROPIC CRITICAL EXPONENTS 

We shall perform the computation of the critical exponents using the normalization con- 
ditions (previously introduced in Sec. II) for the anisotropic cases. 

We start by using eqs.(2),(3) and (5) in order to determine the normalization functions 
Z^( T ), Z^2( r ) in powers of u T along with the Feynman integrals at the loop order required. We 
employ the orthogonal approximation in the calculation of Feynman integrals. The resulting 
expressions for the normalization functions can be written in terms of the four diagrams 
72,-^3,-^4, 15 of one-, two- and three- loop order, respectively. They are presented in Appendix 
A. We recall that each loop integral produces a geometric angular factor ^S m Sd- m T(2 — 
t)^(t)' which shall be absorbed in a redefinition of the coupling constants. Performing 
those redefinitions, we obtain the renormalization functions in terms of the above mentioned 
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diagrams, namely: 

n (^ + 8), 

U 0r = U r [H I 2 U T 

UN + $)I£ f (N* + 6N + 20)I* (5iV + 22)/ 4 , (N + 2)1', 
(iV + 2)/- (iV + 2)(iV + 8)(/ 2 /--4) 2 

%r)= 1 + Ur + ^ ' ( 14b ) 

(iV + 2)I 2 

^0 2 (r) = J- H g «T 

+[ (^ + 7iV + 1 0)/ 2 _^ ( (iV + 2) i | +/4)lii? (Uc) 

As usual, the coefficients of the various terms in powers of u T have poles in e^. These poles 
cancel in the calculation of the (3 T and the critical exponents. In fact using Eqs.(5), the (3 T 
and Wilson functions in either subspace are given by: 

Pt = —T€ L u T [l - a lT u T + 2(a\ T - a 2T )u 2 T ], (15a) 
7<Kr) = -re L M T [26 2 rMr + (3b 3 r - 26 2 rai r )^], (15b) 
70 2 (r) = re L u T [c lT + (2c 2r - c\ T - a lT c lr )u T ]. (15c) 

Substitution of the explicit values of the coefficients given in (13) along with the results for 
the diagrams computed in Appendix A yields the following expression for (3 r 

[3 T = -TU T [e L - — ^— "(! + (Mm - l))«r 
(3JV + 14) , 4 
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<]+0«). (16) 



From this expression, a zero of order of each (3 T characterizing the independent non- 
competing (r = 1) and competing (r = 2) subspaces correspond to the same value of the 
coupling constant (itioo = u 2oo = Uoo), i-e., 



6 

e L <l + e L 



8 + N 



^ (9iV + 42) 



(8 + A^ 



(17) 



Replacing this value back in the functions 7^( T ) and 7</>2( T ) together with the coefficients given 
in (14) and in Appendix A, we find 

/ % r 2 N + 2 . ,6(3AT+14) 1.. 
7«r)(«oo) = ^ T^pt 1 + e ^ (A r + 8)2 - 4 )] " (18) 
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Notice that these expression are exactly the same of those for the critical exponents r]L2 = 
7/1 and 7?£ 4 = r]2 obtained previously in [10j upon the identification 7^(1) (woe) — VL2 and 

7^(2) («oo) = %4- 

The exponents z/ r can be found analogously from the expression of the anomalous di- 
mension of the composite operator 4> 2 . Recalling that 7^,2 t T \ = 7^2 m + 7^( r ), it follows that 
v~ x = — d^2( r ) = 2r — 7<£2( T )(ii TOO ) — 7^( r )(«oo)- Using the results in Appendix A once again, 
we obtain the following intermediate step: 

702(t)(m t ) = T ^ N q 2 ^ t[1 + ei([z 2 ] m - 1) - ^«t]. (19) 

Then, inserting the value into the last expression as well as 7^,( T )(w 00 ) from Eq.(18), we 
get 

1 (N + 2) e 1 (iV + 2)(A 2 + 23A + 60) 2 
" T " 27 + 4r(iV + 8) eL + 87 (ivTs? ^ (20) 

These exponents are in perfect agreement with those calculated previously using a massless 

theory [lCj when we identify v\ = u^ 2 and u 2 = ^L4, respectively. The other exponents can 

be obtained from these (computed diagrammatically) via scaling relations. 



V. ISOTROPIC CRITICAL EXPONENTS 



We are going to use the main results derived in the previous section for the normalization 
functions in the presence of a single mass parameter and only one type of length/momentum 
scale corresponding to the isotropic critical behavior. In this case, each loop integral produces 
an angular factor of S m , the area of an m-dimensional unit sphere, which shall be absorbed in 
a redefinition of the coupling constant. First, we describe the exponents using the orthogonal 
approximation for Feynman diagrams. Then, we shall present the results without the need 
for approximations in calculating the graphs, which is a simple feature of isotropic points. 



A. Exponents in the orthogonal approximation 

The basic definitions of the dimensionless renormalized coupling constant and normal- 
ization functions can be retrieved from Eqs.(14) using the subscript r = 3 in order not 
to confuse with the competing sector of the anisotropic cases. Moreover, the Feynman in- 
tegrals appropriate for the isotropic case have to be used. The results of these diagrams 
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using the orthogonal approximation are given in Appendix B, which are going to be used in 
this subsection. Furthermore, the functions obtained from the normalization constants and 
dimensionless coupling constant through Eqs.(6) in the isotropic cases now read 

As = -€lUs[1 - a 13 u 3 + 2(a 2 3 - a 23 )ul], (21a) 
7^(3) = -e L u 3 [2b 23 u 3 + (36 33 - 2623^13)^3], (21b) 

70 2 (3) = £L%[Cl3 + (2C 23 - cj 3 - Oi3Ci 3 )tt 3 ]. (21c) 

Recall that in the above expressions = 8 — m. Using the integrals presented in Appendix 
B, we find 

(N + 8) n 1 
P3 = -u 3 [e L (1 - -eL)u 3 

+ M+M f41 + 0(t a. (22) 

From the eigenvalue condition ,53(1*300) = we obtain the following result for the zero of 
order of the dimensionless coupling constant 



6 J 1 M tL 



1 (9N + 42) 
2 + (8 + iV) 2 



(23) 



Using this value in the expression for 7^(3) in conjumination with Eqs.(14) and the results 
in Appendix B, one can show that 

2 N + 2 ^ ,6(3iV+14) 1 N1 
7«*)(«aoo) = e L 4(jV + 8)2 [l + e L ( (iV + g)2 - -)]. (24) 

Once more, this is precisely the value of the exponent t\t A previously obtained using the 
renormalization group equation in the massless theory [lO]. Moreover, using again the 
results in Appendix B in the definition of 7<?i 2 (3); we obtain 

(N + 2) 1 1 

7^(3)(«3j = g u 3 [l - -e L - -u 3 \. (25) 

When Eq.(23) is replaced into this expression it yields 

x (N + 2) r (3iV + 9) . 
7^(3)(%oo) = 7^vT8) eL[1 + (N + 8)* €l] ' (26) 

and recalling that v 3 = v L4 = (4 - 702 { 3 ) (w 3oo ) - 7<H3)(«3oo))~ 1 , we get to 

_ 1 (N + 2) 1 (AT + 2)(iV 2 + 23iV + 60) 2 

^ L4 "4 + 16(iV + 8) eL+ 256 (7V + 8)3 ^ ^ 
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The expression for 70 2 (3)(' u 3oo) in Eq.(26) is the same as the one associated to a scalar theory 
in the massless limit, computed at the fixed point using the renormalization group equation. 
Besides, Eq.(27) corresponds to Eq.(208) for this exponent in the orthogonal approximation 
using the method of Ref . jlo| . 



B. Exponents in the exact computation 



The Feynman integrals of the isotropic case can be manipulated without referring to any 
approximation in the massive case, as described in Appendix C. Our goal here is to rederive 
the exponents using the machinery developed in the last subsection replacing, however, the 
expression for the diagrams by their exact counterparts presented in Appendix C. 

We start by writing down the (3$ function explicitly in terms of the diagrams calculated 
in Appendix C, which in that case reduces to 

(N + 8). 1 , 



-u 3 \e L 



6 



4 



(41A + 202) 2 4 
-u 3 \ + 0{u 3 ) 



1080 



(28) 



Therefore, this expression has a zero at the following value of the dimensionless coupling 
constant 

6 (\ e L \l (A1N + 202)1 ] 

{ ZJ > 



U 3 c 



s 1 N 6L V +£ i 



(41N + 202) 



15(8 + iV) 2 

Following the same steps as before it is not difficult to show that this value of the coupling 
constant yields 



7<K3)(«3c 



— e, 



3(iV + 2) 



+ 



(N + 2) r (41iV + 202) 23, 



(30) 



:L 20(iV + 8) 2 ' 10(iV + 8) 2L 10(iV + 8) 2 ' 80 J 

This expression is equal to the value of the exponent T]L4 calculated exactly in conformity with 
that of 0{e\) firsty found by Hornreich, Luban and Shtrikman \X\ and recently confirmed 
and extended to 0(e|) in Ref. Jj]] for the massless case, as expected. 

Finally, let us determine ula- First, we have to make use of the results in Appendix C to 
write down the following expression for arbitrary dimensionless coupling constant 



- f \ (^ + 2 ) n 1 1 1 

7^(3) («3j = g Ms[l - + ^U 3 ] 



(31) 
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Second, we replace the especial value of the coupling constant u 3oo in this equation in order 
to show that 

/ \ (N + 2) (13JV + 41) . , oox 

7^(3) (tiacc) = ^Ts)^ 1 - i5(iV + 8)^ i] ' (32) 

which allows us to figure out the value of vn from 7^3) (u 3oo ) and 7^2(3) (w 3oo ), namely 

1 (N + 2) (A + 2)(15A 2 + 89iV + 4) 2 

UL4 ~ 4 + 16(iV + 8) eL + 960(iV + 8)3 ^ {66) 

I — I 

This is precisely the exponent in the exact calculation recently found at 0(e\ [lSj. As 
before, the other exponents can be discovered through the use of scaling relations. 



VI. DISCUSSION OF THE RESULTS AND CONCLUSIONS 

It is interesting to find two different methods producing the same values of the criti- 
cal exponents calculated diagrammatically. Those previously computed using a massless 
scalar field theory subtracted at nonzero external momenta using the renormalization group 



equations in Refs. 10| are rederived herein using massive theories (with two or one inde- 
pendent mass scales depending whether we treat the anisotropic or isotropic cases) in the 
Callan-Symanzik method. 

n 

Most of the remarks done originally by Brezin, Le Guillou and Zinn- Justin in Ref . |14| 
for conventional critical phenomena are still valid here, for isotropic as well as anisotropic 
cases. First, the eigenvalue condition (3 T {u TOQ ) = corresponds to a scale invariant theory, as 
shown here, but it is not attractive: the theory is only scale invariant precisely at this point, 
since the derivative j3' T {u roo ) calculated at this point is positive. Second, our calculations are 
valid in the region of momenta much larger than the mass in any particular subspace. The 
difference dwells in the anisotropic cases: two inequivalent subspaces require two independent 
mass scales along with two coupling constants which consistently have the same eigenvalue 
leading to scale invariance separately in each subspace. Universality is recovered when the 
renormalized dimensionless coupling constants are fixed at the solution of the eigenvalue 
u too . 

One step further within the formalism introduced in the present work is to extend it 
in order to treat arbitrary competing systems, when arbitrary types of competing direc- 
tions characterized by higher order derivative terms of arbitrary powers are involved [is)], 
generalizing the present concept. The Feynman integrals get much more difficult to solve, 

15 



especially when one tries to perform the exact calculation for the isotropic cases. That is 
the reason we have calculated here the integrals using the orthogonal approximation and 
the exact computation for those integrals. In spite of yielding different values for the critical 
exponents, the difference among them in each calculation is negligible as already pointed 
o„U„Q. We shalll e a ve th e subject of generic competing systems in the Callan-Symanzik 
method for future work 

The anisotropic behavior with two independent mass scales and the analysis developed 
herein might be useful in addressing further anisotropic problems in quantum field theory, 
possibly with broken Lorentz/Poincare invariance. First, it is useful to consider the higher 
order terms in momentum space appearing isotropically in the ultraviolet regime. This pos- 
sibility arises, for instance, in Lorentz-violating dispersion relations for cosmic ray processes 
due to the nontrivial character of the short-distance structure of spacetime 17|. This effect 



can be implemented in the propagation of massive particles as well, even anisotropically. 
One mathematical possibility to describe such phenomena is that the metric tensor depends 
upon the momentum variables. This has been done recently to investigate the simplest ex- 
ample in anisotropic Planck-scale modified dispersion relations for a relativistic particle of 
mass m propagating in 2 dimensions taking into account phenomenological quantum gravi- 
tational effects (i.e, the modified dispersion relations) in the framework of Finsler geometry. 
The coefficient of the cubic term in the momentum pi in the mass-shell condition is directly 
related to the Planck mass M [18]. In principle, these modified dispersion relations can be 
postulated for massive fields as well. Then the procedure should be the opposite to that de- 
scribed in the present paper to obtain the Lagrangian: keep So 7^ 0, restore the dimensionful 
character of the higher order derivative terms and introduce the equations of motion with 
any power of momenta (mass-shell condition) isotropically or anisotropically, as a constraint 
using an auxiliary field (Lagrange multiplier). The problem is that the equation of motion 
obtained by varying the Lagrange multiplier is not easy to solve (even in the particle case) 
and the renormalization group structure in this framework is not obvious. 

Another perspective is to follow the trend proposed here of maintaining 5q = which 
allows the decoupling of the inequivalent subspaces and makes it easier the ultraviolet renor- 



malization of quantum fields. Sie, 
Poincare invariant field theories 



gel's recipe to introduce masses by dimensional reduction in 
19| would be modified in the present case with at least two 



inequivalent subspaces. This would require two extra dimensions to characterize completely 
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the two mass scales involved in the problem, one for each momentum subspace. This is so 
because one can mathematically formulate a gradient operator for "Lifshitz spaces" which 
can be defined as <9; = (<9i, 9newW=m+i) = 9 L(d-m+i)' -i d new(d) = d 2 old{d) ) through 
the dimensional redefinition proposed in lOj, with the same Euclidean metric, keeping in 



mind that the mass mi characterizes the noncompeting (quadratic) subspace, while m 2 is 
the typical scale of the competing (quartic) subspace, with the latter having half of the 
canonical dimension of the former. This new operator would reflect directly the "folding" of 
the dimensions along the competing subspace as a dynamical effect provoked by the combi- 
nation of attractive and repulsive components of the background on which quantum fields 
propagate. This new flat space limit does not correspond to Minkowski space. Although 
this is trivial for scalar fields, it is definitely not the case for defining spinors and vector 
fields in these Lifshitz spaces with anisotropic higher derivative terms in the kinetic part. 

In conclusion, the new method proposed in the present work to study m-axial Lifshitz 
points for the massive scalar field theory renormalized at zero external momenta opens the 
possibility to study field theories with inequivalent subspaces characterized by two (or more) 
mass scales. Universality is confirmed, since the critical exponents calculated in the massive 
theory using the Callan-Symanzik formalism agree with those obtained in the massless scalar 
field theory renormalized at nonvanishing external momenta. 
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APPENDIX A: FEYNMAN GRAPHS FOR ANISOTROPIC BEHAVIORS 

We compute the appropriate diagrams corresponding to one-, two- and three- loop inte- 
grals using dimensional regularization and Feynman parameters to fold several denominators 
together. Only the one-loop integral for the four-point function can be solved exactly for 
anisotropic criticalities. The higher loop diagrams can be calculated, however, using the 
orthogonal approximation, which is the most general approximation consistent with the 
physical property of homogeneity. The basic statement of this approximation is that the 
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loop momenta in a given subdiagram is orthogonal to all loop momenta appearing in other 
subdiagrams. Due to our normalization conditions defined in the text, we shall work with 
zero external momenta. Explicitly, the integrals we have to resolve are given by: 



= r d d - m gd m k 



I' 3 = dh Qp2 K) \p=K=o(= ^£P-\p=k=o) where h(P,K) is the integral 

d d - m q 1 d d - m q 2 d m k 1 d m k 2 
(ql + (kfy + 1) {ql + (P 2 y + 1) [( qi + q 2 + Pf + {(h + k 2 + Kf) 2 + 1] ' 

d d - m q 1 d d - m q 2 d m k 1 d m k 2 



r d d - m q x d d - 
J (q\ + (k\f + 1 



f(q 2 2 + (klf + l) 



-L 

-1 , , 7, I 

77 



where [i 2 \ m = 1 + \(^(l) - ij>{2 - f )) and 

ln[x(l — x)] 



C 1 

7o 1 - ^(i - a:) + 



(A2) 



x 5 , (A3) 

\(qi + qrf + {(h + k 2 Y)\l] 

and = 9Ib Qp2 K ^ \p=k=o(— 9Ir °dK^ \p=k=o) where I$(P, K) is the three-loop diagram 

d d - m qi d d - m q 2 d d - m qzd m kid m k 2 d m h 
(( Ql + Pf + ((h + Kff + 1) (ql + (MP + 1) (ql + (fc 3 2 ) 2 + 1) 

x 5 3 . (A4) 

[(«i + ? 2 ) 2 + ((*i + ^ 2 ) 2 ) + l][(?i + qs) 2 + ((fci + ^ 3 ) 2 ) + 1] 

Proceeding as indicated above, their e^-expansion have the following expressions 

J 2 = l[l + ([2 2 ] m -l)e L ] + 0(e L ), (A5) 

I's = ^[1 + (2[<2]m " " ^ + O(cl), (A6) 

^=2^-(l + (2[<2]m-^L + 0(4)), (A7) 



Js = 64 [1 + (3Nm " 4 )£L + ° (e ' ] " 4^ 7 ' (A8) 



[1 -x(l -x)f" (A9) 
is the same integral appearing in the original work by BLZ. This integral disappears in 

the calculation of the critical exponents. Its appearance here is a general feature of the 

orthogonal approximation, since it will be present in the final result of the calculation of the 

diagrams in the isotropic cases as well. It does not appear in the masless theory and is only 

an artifact of our choice of the normalization point for the external momenta. 
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APPENDIX B: ISOTROPIC DIAGRAMS IN THE GENERALIZED ORTHOG- 
ONAL APPROXIMATION 



We shall now turn our attention to the isotropic integrals in an expansion in e L = 8 — d. 

The simplifying feature of this case is the existence of only one type of correlation length, 

which implies only one mass as usual in field theory calculations. The issue here is a quartic 

power in the propagator. The expressions for these graphs are: 

f d m k . > 

h = / 5 , (Bl) 



I' s = 9l3 Q^'I < ^ \k=o where h(K) is the integral 

d m k 1 d m k 2 



{{k\f + 1) {{k^ + 1) [({h + k 2 + Kf) 2 + 1] ' 
d m k 1 d m k 2 



(B2) 



h = 



my + if my + 1) 



l({k 1 + k 2 yy+i] 

and I' 5 = dI Q^V \k=o where h{K) is the three-loop diagram 



x— * 2 ~ , (B3) 



d m k x d m k 2 d m k z 



+ Kyy + 1) {{kjy + 1) {{kjy + 1) 

x 5 5 • (B4) 

[((^i + ^) 2 ) 2 + i][{{h + k 3 y) 2 + i] 

Using the same steps as indicated in the orthogonal approximation for the anisotropic be- 
haviors we find for their e^-expansion the following results 

h = -[l-je L ]+0{e L ), (B5) 

t-^p-r-i-k'+o^ (B6) 

/ 4 =2^(l-j£L + 0(4)), (B7) 

^"Ml 1 -?^ ^-^ (B8) 
where J is the same integral which appears in Appendix A. As before this integral do not 

contribute for the critical exponents. 
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APPENDIX C: ISOTROPIC INTEGRALS IN THE EXACT CALCULATION 



We shall consider now the integrals without performing approximations. The expressions 
are the same of those given in the previous appendix. Before quoting the results some 
commentaries are in order. First, the simplicity of the isotropic cases allow us to perform 
the integrals exactly as advertised. A useful trick is to write the propagators in the form 
p2js+ij — (k 2 +i)(k 2 -i) an d then use a Feynman parameter. One can proceed as before 
since now the propagators are quadratic in the loop momenta, but the resulting parametric 
integrals turns out to be more complicated than those in the orthogonal approximation. 
Moreover, some care must be exercised in order to manipulate the remaining factors of 
i. Setting = in the powers of i when performing the parametric integrals effectively 
suppress it. We can then figure out the final expressions for the exact calculation, namely 

h = -[l-^ L ]+0(e L ), (CI) 

'^^-w^-\ H+0 ^ < C2 > 

^ = ^-(i-^L + 0( e 'i)), (C3) 

^^-i^ ^!-^ < C4 > 

and now the multiparametric integral is specific to m-axial isotropic Lifshitz points and 
turns out to be 

-1 r l /•! /•! rl 



-V77 1 

Jo Jo Jo 



H = / dzz(\ — z) I I I dxdydt I duu{\ — u 



\3 



x ln(u [ 2 " (1 ~ " } + 2 f" ~ 1 - 2t + 1] + 2t - l) . (C5) 
V z[l — z) J 

The integral H (as well as the integral I appearing in the orthogonal approximation) does 
not contribute to the critical exponents since it disappears in the renormalization algorithm. 
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